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Abstract: Soil salinization is a major problem affecting soils and threatening agricultural sustainability in 
arid and semi-arid regions, which makes it necessary to establish an efficient strategy to manage soil salinity 
and confront economic challenges that arise from it. Saline soil recovery involving drainage of shallow 
saline groundwater and the removal of soil salts by natural rainfall or by irrigation are good strategies for 
the reclamation of salty soil. To develop suitable management strategies for salty soil reclamation, it is 
essential to improve soil salinity assessment process/mechanism and to adopt new approaches and 
techniques. This study mapped a recovered area of 7200 m?to assess and verify variations in soil salinity in 
space and time in Kairouan region in Central Tunisia, taking into account the thickness of soil materials. 
Two electromagnetic conductivity meters (EM38 and EM31) were used to measure the electrical 
conductivity of saturated soil-paste extract (ECe) and apparent electrical conductivity (ECa). Multiple 
linear regression was established between ECe and ECa, and it was revealed that ECa-gm3sis optimal for 
ECe prediction in the surface soils. Salinity maps demonstrated that the spatial structure of soil salinity in 
the region of interest was relatively unchanged but varied temporally. Variation in salinity at the soil surface 
was greater than that at a depth. These findings can not only be used to track soil salinity variations and 
their significance in the field but also help to identify the spatial and temporal features of soil salinity, thus 
improving the efficiency of soil management. 
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1 Introduction 


Soil salinization is one of the major reasons for agricultural land degradation, especially in arid 
and semi-arid regions (Marlet, 2004). In arid and semi-arid regions across the world, 
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irrigation-dependent agriculture is a threat to the environment as it can cause soil degradation by 
salinization or sodization (Huang et al., 2015). Pitman and Lauchli (2004) reported soil 
salinization to be widespread in irrigated agriculture, particularly in Asia and Africa, where a 
heavily affected country is Tunisia. 

In Tunisia, salinization of agricultural land constitutes a major problem as the area affected by 
salinization covers approximately 1.5x10° hm’, which is about 10% of the total area of the 
country. The affected area is located mainly in the center and south of Tunisia. Soil salinity maps 
can aid in better evaluation and monitoring of the extent of salt accumulation in any given area, 
and highlight regions requiring reclamation. To address land degradation caused by salinization, 
there have been attempts to restore and develop several thousand hectares of saline and 
hydromorphic soils (Hachicha, 2007). 

To study salt-affected soils, researchers have developed several approaches and techniques to 
measure the spatial and temporal variability in the salinity of soils and the natural and 
anthropogenic factors of dependence associated with them (Corwin and Lesch, 2003; Hachicha, 
2007). Since the early 1950s, the United States Salinity Laboratory Staff (USSL) has standardized 
salinity assessment and its interpretation. Originally, electrical conductivity (EC) was used to 
determine electrical conductivity of saturated soil-paste extract (ECe) or other aqueous extracts of 
soils. However, these methods have proved insufficient. Furthermore, they are expensive, require 
time and effort, and can only be performed on a limited number of samples, which limit their 
usefulness for the spatial and temporal monitoring of soil salinity (Corwin and Lesch, 2003). 
Alternative measurement methods have been developed since the 1960s, one of which is the 
measurement of apparent electrical conductivity (ECa). ECa instruments have been widely used 
to map the spatial distribution of soil salinity in Canada and Australia (Thomas et al., 2009). 
These instruments have been increasingly used in the agricultural development of countries, such 
as Senegal (Job et al., 1995; Barbiero et al., 2001), Morocco, Spain, Tunisia, Turkey (Aragiiés et 
al., 2011), Niger (Adam et al., 2012), and Algeria (Berkal et al., 2014). Tools to measure ECa are 
used to not only map the spatial distribution of salinity (Corwin and Lesch, 2003) but also 
characterize soil salinity at multiple depths (Paz et al., 2020; Sousa et al., 2021; Wang et al., 
2021). There are many commonly used instruments. The best known are the EM31 (EM, 
electromagnetic conductivity meter), EM34, and EM38 marketed by Geonics Ltd. (Mississauga, 
Canada), DUALEM-1 (Dualem Inc., Milton), and DUALEM-2 (Dualem Inc., Milton). 

In salt-affected soils, salt dominates the response of the EM and, generally, there is a good 
correlation between ECa and ECe (Rhoades et al., 1989; Amakor et al., 2014; Cassel et al., 2015). 
Relative values obtained from ECa measurements can be calibrated against known EC values in 
the soil solution and EC measurements in the laboratory (Rhoades et al., 1989; Bellier et al., 1997). 
This calibration may lead to biases due to in situ sampling and differences in water content and 
texture of the soil volume used (Johnson et al., 2005). To characterize the spatial and temporal 
variability in soil properties (mainly salinity), researchers developed a geophysical method 
combining geostatistics and point measurements of EM induction (Gascuel-Odoux and Boivin, 
1994; Job et al., 1995; Hendrickx et al., 2002; Michot et al., 2003; Wang et al., 2021) that allows 
quick, non-destructive, and precise estimations of salinity at different depths (Triantafilis et al., 
2001a; Padhi and Misra, 2011). Numerous studies on soil salinity have shown that spatial estimates 
based on single-point measurements of EC have limited accuracy (Hajrasuliha et al., 1991; 
Hosseini et al., 1994; Odeh et al., 1998; Walter et al., 2001; Douaoui et al., 2006) due to soil 
sampling may not be adapted to account for small-scale soil salinity changes as well as 
non-stationary salinization processes. Technical capabilities to study soil salinity have now been 
expanded, by remote sensing, geographic information system (GIS) using geostatistical methods, 
digital soil mapping approaches, simulations conducted by modeling software, and many other 
modern instrumental methods are available for assessment and mapping of soil salinity. 

The most widely used approach is the establishment of linear regression calibrations, which can 
predict the average ECe at different depths (e.g., <1.2 m) and ECa (Amezketa and de Lersundi, 
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2008). However, this approach leads to information loss regarding salt content at specific depths 
as EM readings were performed in two modes, i.e., two investigation depths, depending on the 
devices used (EM38 or EM31). Williams and Baker (1982) highlighted EM that can access data 
on salinity from greater depths, and thus better delineate areas that may become saline upon salt 
mobilization by increased groundwater levels and land-use changes. 

Many soil salinity studies have explored the regional distribution of saline soils and the 
characteristics of salinization, yet few data are available on vertical EC distribution and the 
transfer of salt to deeper soil layers. In semi-arid regions, saline soil is remediated by water 
ponding or irrigation, which improves soil fertility (Provin and Pitt, 2001); natural desalinization 
may occur due to the occasional intensive rainfall that transfers large amounts of salts to deeper 
soil layers, increasing the salinity of groundwater (Trabelsi et al., 2005). It is thus crucial to 
include deep soil layers in the assessment of soil salinity. Therefore, our study aimed to couple 
two EM devices, EM31 and EM38, to monitor and describe the spatial and temporal variations of 
soil salinity in deeper soil layers using multiple linear regression (MLR) models on a study site 
that was being used as an orchard with irrigated cultivation of pomegranates plants in Kairouan 
region, Central Tunisia. 


2 Materials and methods 


2.1 Study area 


The study area (35°48'33"N, 10°02'35"E) is located about 15.0 km north of the city of Kairouan 
in Central Tunisia. It is approximately 4.5 km long and 1.0 km wide and covers an area of about 
450 hm?. It is geographically located in the alluvial plain of Bas Sisseb-El Alem and 
administratively situated in the Governorate of Kairouan and the Delegation Sbikha. 

The region of interest (represented by the transect in Fig. la) comprises the alluvial plains of 
Sisseb-El Alem that are characterized by contrasting morphology of the borders of the plain, where 
the wind and water erosion processes predominate, and eroded material from the plain accumulates 
in particle size (Safar, 1983). The soils are generally stratified, and the alluvial clay at the surface is 
mostly dominated by montmorillonite mixed with kaolinite, illite, and traces of quartz (Safar, 1983). 
This is an upper arid climatic zone that is strongly influenced by winds coming from the 
Mediterranean Sea. The average annual precipitation is about 250.0 mm, with 120.0 mm falling in 
autumn, 40.0 mm in winter, and 63.0 mm in spring. The annual mean temperature is 17.9°C, with a 
maximum monthly mean of 37.8°C in July and a minimum monthly mean of 10.6°C in January. 
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Fig.1 An aerial view of the study area (a) and location of soil sampling profiles (b) 
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The second region of interest (prospected site in Fig. 1b), is about 90.0 m long and 80.0 m wide 
(7200 m°) and is located at 35°48’32"N and 10°02'35"E. The region is cultivated with 3-year-old 
pomegranate plants in 20 rows with an inter-row distance of 4 m and a space of 2 m between each 
plant. 


2.2 Data collection, sampling and analysis 


The EM31 and EM38 use EM technology to measure ECa in soils without planting electrodes. 
They consist of two coils, a transmitter, and a receiver coil that are spaced 3.7 m apart in EM31 
and 1.0 m apart in EM38. The device can be used in two vertical dipole modes in fact the 
theoretical survey reaches 1.2 m by EM38 and 6.0 m by EM31 and also two horizontal dipole 
modes whose theoretical survey of order 0.8 m for the EM38 and 3.0 m for the EM31. Three 
campaigns to collect ECa measurements were carried out in spring (April) 2015, autumn 
(October) 2015, and autumn (October) 2016 using both meters. ECa values were then collected 
from the spacing between 20 rows of pomegranates, where measurements were made so that a 
single measurement of four plants could be obtained. 

A transect (section) leads from the entrance of the study area to the west toward the river bed of 
Oued Bogal, which is present about 4.5 km to the east of the study area. Calibration profiles were 
selected randomly. 

Each sampling region or transect consisted of 218 single values collected as follows. During 
each geophysical measurement, 23 sites (profiles) were sampled with the EM38: 15 sites 
(profiles) were sampled from prospected site to analyze the spatial variability in salinity within 
the plot and eight sites were sampled along the 4.5 km transect to measure regional variability in 
soil salinity. EM31 measurements were collected at six sites: three sites were within prospected 
site and three sites were located along the transect (Fig. 1). 

At each EM38 measurement site, soil samples spanned in the range of 0.2—1.2 m from the top. 
At the EM31 measurement sites, soil samples spanned in the range of 0.2 to 6.0 m from the top. 
For each sample, gravimetric water content (GWC) and ECe were measured as described 
previously by Rhoades et al. (1989). The average ECe for EM38 measurement sites was 
calculated for two soil layers: 0.0-0.6 m (surface soil layer) and 0.6—1.2 m (subsurface soil layer). 
The EM31 measurement sites were divided into four layers: 0.0—0.6 m (surface soil layer), 
0.6-1.2 m (subsurface soil layer), 1.2-3.0 m (intermediate deep soil layer), and 3.0—6.0 m (deep 
soil layer). 

Four ECa measurements were recorded. Specifically, two measurements were recorded with 
the EM38, where the first was in the horizontal dipole mode (ECa-em3sn), i.e., measured with the 
EM38 coils positioned vertically for an investigation depth of 0.8 m and the second in the vertical 
dipole mode (ECa-em3sy), i.e., with the coils positioned horizontally for an investigation depth of 
1.2 m (McNeill, 1980). The other two measurements were recorded using an EM31 in vertical 
(ECa-Em3iv) and horizontal (ECa-gm3in) dipole modes for penetration depths of approximately 6.0 
and 3.0 m, respectively (McNeill, 1980), with the instrument held 1.0 m above the soil surface. 
The EM31 had an inter-coil spacing of 3.7 m and was operated at a frequency of 6.4 kHz. Each 
ECa-em (apparent electrical conductivity measured by the EM (dS/m)) measurement location was 
georeferenced using the 4600 LS™ GPS (with differential correction), which is manufactured by 
the Trimble Italia Srl, Italy. The GPS receiver had a horizontal precision of >2.0 m. 

To compare ECa data (maps) obtained on different dates from soils of varying wetness, 
ECa-emv (apparent EC measured in vertical dipole orientation (dS/m)) and ECa-gmu (apparent EC 
measured in horizontal dipole orientation (dS/m)) values measured with both tools (EM38 and 
EM31) were adjusted to the GWC of 20.00% using the equations of Job et al. (1995) established 
specifically for Tunisian saline silty clay soils: 


ECa-yqui(20%) = BC8-pyagca,) + 7-7 x (20-4), a) 
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EC-pyiy(o0% =EC8-pmvcq) + 7.7 x (20-4), (2) 


where ECa-gmx(20%) (dS/m) and ECa-gmvi20%) (dS/m) are the values of ECa-emn and ECa-gemy 
corrected at a GWC of 20.00%, respectively; @ (%) is the GWC; and ECa-emuca) (dS/m) and 
ECa-emvg) (dS/m) are the values of ECa-emu and ECa-gmy measured at a GWC of 0, 


respectively. These equations can be used to standardize ECa-gm values measured at GWC from 

8% to 29% (Job et al., 1995). 

2.3 Soil salinity mapping 

The extrapolation of salinity values to different depths from the measured sites was carried out 

with statistical models using the EM38 and EM31 datasets. The averages of ECe were calculated 

for the soil layers mentioned before. The objective of this mapping was to estimate the average 

ECe for all non-sampled sites (so) in the study area. The procedure was divided into three steps. 
First, for each calibration, we constructed an MLR model to predict the average ECe for the 

different soil layers according to the method of Rhoades and Corwin (1981): 


ECe = a x ECa-pyy + bx ECa-pyyy +C, (3) 


where ECe (dS/m) is the average EC of saturated soil-paste extract for the soil layer being 
considered; ECa-emv (dS/m) is the apparent EC measured in vertical dipole orientation; ECa-zmu 
(dS/m) is the apparent EC measured in horizontal dipole orientation; and a, b, and c are the 
coefficients of the linear equation. 

Second, the MLR models from the first step were applied to predict the average ECe for 
different soil layers. 

Third, Ordinary Kriging method was applied using the Geostatistical Analyst module of 
ArcGIS 10.5 software (ESRI, 2010) to estimate ECe at all four depths: (0.0-0.6 m (surface soil 
layer), 0.6-1.2 m (subsurface soil layer), 1.2-3.0 m (intermediate deep soil layer), and 3.0-6.0 m 
(deep soil layer)) using p=218 (where p is the number for measurements) predictions of ECeřmLR 


(electrical conductivity predicted by the MLR) from the second step. ECe (sp, h;) was estimated 
using Equation 4 (Webster and Oliver, 1992): 


ECe”(sọ,h,) =J à ECe’ mir (s ;,h,), (4) 

i=l 
where (A, is the weight of the MLR estimate at the i™ site for i=1 to p: neighbors considered in 
the Ordinary Kriging procedure/process; and (sj, hj) is the position of ECe measured. For 


cross-validation, estimates of ECe were compared to the measured average ECe values for each 
soil layer and at n ECe sampling sites by calculating the mean average error (MAE) and 
root-mean-square error (RMSE) (Oliver and Webster, 2014). 

The one-dimensional inversion program, UMR METIS (the METIS Joint Research Unit: 
Environments, Transfers and Interactions in Hydrosystems and Soils), was applied to the large 
mesh (5.0 m) EM data acquired in the site using iterative least-squares inversion method (Guérin et 
al., 1996). Four different datasets were used for each date corresponding to ECa-emssą at 0.6 m, 
ECa-emaev at 1.2 m, ECa-gm3inat 4.0 m, and ECa-gm3iv at 6.0 m. All methods were tested on 218 
monitoring sites at 0.3 m elevation, and data were acquired continuously for two successive 
surveys. 


2.4 Mapping of ECa changes over time 


To evaluate the stability of ECa-rm over time, we calculated the coefficient of variation (CV) of 
measured ECa at each geophysical measurement. 

According to Blackmore (2000), the CV has helped to assess the temporal stability of crop 
yields (Shi et al., 2005), evaluate the changes in soil's chemical and physical properties in 
grasslands (Guo et al., 2015), and examine the variations in soil salinity over time in China. 

The CVs were calculated using EM38 and EM31 datasets measured in both vertical and 
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horizontal modes. The CV of ECa was calculated using Equation 5. 


2 
n9 (ECa; j -— (Eca) 
t=l t=l 


CV; = Z ae , (5) 
[Sec 
t=1 
n 


where CV; is the CV at site i and n is the number of ECa measurements. 
Maps of the CV of ECa were drawn by Ordinary Kriging using the Geostatistical Analyst 
module of ArcGIS 10.5 software (ESRI, 2010). 


3 Results 


3.1 Soil GWC and soil salinity profiles 


The variation in GWC is illustrated in Figure 2. The GWC profiles had a similar shape for the 
three sampling periods (Fig. 2). The GWC varied from 21% to 44% in the surface soil layer 
(0.0-1.2 m). Although salinity levels were very high in most ECe profiles, they varied with time. 
The ECe ranged from 11.18 to 29.80 dS/m in the surface soil layer with an average value of 13.79 
dS/m, while GWC varied from 9% to 29% in the deeper soil layers. The soil surface was wetter in 
autumn 2016. An ECe of 66.50 dS/m was recorded below 3.0 m from the surface, and 40.0 dS/m 
from 4.0 to 6.0 m. Slight surface desalinization was observed from spring 2015 to autumn 2016, 
while salinity increased with depth in all seasons. 


3.2 Correlation between ECa-gmvy and ECa-rmn measurements 


The results of regression analysis between ECa-gmvy and ECa-emu for all calibration sites are 
shown in Figure 3. There was a clear linear relationship between ECa-emy and ECa-gmu at all 
sampling sites (29 sites) using EM31 and EM38 datasets. All ECa-gmy values were higher than 
ECa-emu values, indicating an increase in soil salinity with soil depth. 


3.3 Seasonal variation in ECa-gm measurements 


Table 1 presents descriptive statistics of EM31 and EM38 measurements for each sampling. 
ECa-em measured by the EM38 was always lower near the surface than at a depth 
(ECa-Em3gv/ECa-Em3sn> 1.00). Mean ECa-gm3sy values were higher during spring 2015 (6.97 
dS/m) and decreased significantly in autumn 2015 (6.49 dS/m), which is attributed to rainfall that 
preceded the measurements. Finally, the mean ECa-gm3gy returned to its original value in autumn 
2016 (6.97 dS/m). The ECa-gm3gv measurements showed the highest significant variability in 
autumn 2015, with a minimum value of 3.71 dS/m and a maximum value of 10.11 dS/m. The 
mean ECa-em3su value was 4.31 dS/m in spring 2015, which decreased progressively to 3.88 dS/m 
in autumn 2015 and to 3.42 dS/m in autumn 2016. The CV of ECa-zms3gv decreased from 15.67% 
in spring 2015 to 11.36% in autumn 2016. During the monitoring period, the lowest CV of 
ECa-gm3iv value (4.30 dS/m) was recorded in autumn 2016, while the highest value of 7.95 dS/m 
was recorded in spring 2015. The ECa-gm3su values were characterized by low CV ranging from 
16.70% in spring 2015 to 24.10% in autumn 2016. 

The ECa-ev3iv and ECa-em3in values measured with the EM31 were lower than those 
measured with the EM38. The mean ECa-gmaiy value of 4.54 dS/m was recorded in spring 2015, 
which remained constant until autumn 2015 and then increased to 5.36 dS/m in autumn 2016. A 
minimum value of 3.79 dS/m was observed from spring to autumn 2015 and the highest value of 
5.97 dS/m was measured in autumn 2016. The mean ECa-gmain values were close to 2.09 dS/m 
from spring to autumn 2015; they increased to 3.54 dS/m in autumn 2016 (Table 1). 
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Fig. 2 Seasonal variations in different salinity-related metrics recorded in spring 2015, autumn 2015, and 
autumn 2016. (al—a2), gravimetric water content (GWC) and electrical conductivity of saturated soil-paste 
extract (ECe) in the soil layer of 0.0-6.0 m; (b1—b6), GWC and ECe in the soil layer of 0.0-1.2 m. P1-P5 
represent five profiles, which respectively contain five soil layers: 0.0—0.2, 0.2—-0.4, 0.4-0.6 m, 0.6-0.8 and 
0.8-1.2 m. 
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Fig. 3 Plot of apparent EC measured in vertical dipole orientation (ECa-emv) versus apparent EC measured 
in horizontal dipole orientation (ECa-rmu) at all sampling sites (29 sites) using EM31 and EM38 datasets 


Table 1 Descriptive statistics of ECa-emu and ECa-emv in spring 2015, autumn 2015, and autumn 2016 


Spring 2015 Autumn 2015 Autumn 2016 


Statistic EM38 EM31 EM38 EM31 EM38 EM31 


Eca-emv Eca-emu Eca-emyv Eca-emu Eca-emy Eca-emy Eca-emv Eca-emu Eca-emy Eca-emu Eca-emy Eca-emu 


Mean (dS/m) 6.97 431° 4.54 3.02 6.49" 3.88* 4.52 2.96 6.96 3.42” 5.36* 3.54" 
Med. (dS/m) 6.92 4.24 4.48 2.95 6.46 3.77 4.47 2.93 6.93 3.36 3:37 3.51 
Min. (dS/m) 4.58 2.84 3.79 2.09 3.71 2.63 3.79 2.48 4.42 2.23 4.65 2.91 
Max. (dS/m) 9.35 6.59 5.45 4.46 10.10 5.81 5.45 3.91 9.14 6.86 5.97 4.58 
SD (dS/m) 1.09 0.72 0.36 0.40 1.01 0.90 0.35 0.27 0.79 0.82 0.23 0.24 
Kurtosis —0.86 0.14 —0.30 2.56 0.39 0.18 0.29 0.42 0.38 1.79 1.62 0.16 
Skewness 0.01 0.54 0.34 1.27 0.34 0.02 0.38 0.74 0.13 1.22 0.72 0.02 
CV (%) 15.67 16.70 7.95 13.28 15.50 23.17 7.83 9.19 11.36 24.10 4.29 6.77 


Note: ECa-gmv, apparent EC measured in vertical dipole orientation; ECa-gmu, apparent EC measured in horizontal dipole orientation; 
Med., median; Min., minimum; Max., maximum; SD, standard deviation; CV, coefficient of variation. * indicates a significant difference 
from average values at the 0.05 level according to the Tukey test. 


Moreover, the lowest ECa-em3in value of 2.09 dS/m was measured in spring 2015, while the 
highest value of 4.58 dS/m was measured during the autumn 2016. The ECa-em3iv and ECa-em3in 
data collected during the three sampling periods with the EM31 exhibited a moderate variability 
as indicated by their low CVs of 4.29%-13.28%. 

The variation should be considered stable when CV is less than 10.00%. Also, the CV is stable 
between 10.00%—25.00% and is unstable when CV is higher than 25.00%. In addition, the mean 
ECa-em31v/ECa-gm3in ratio remained about 1.50 during the three sampling periods, indicating that 
soil salinity was higher in the deeper soil layers and there were very few changes in soil salinity 
over time between the different sampling depths. 

Conclusively, the variation in ECa-emssą was significant in the three sampling periods, whereas 
the variation in ECa-em3sy was significant for autumn 2015 only. Further, for the EM31 
measurements, variations were significant only for autumn 2016. 


3.4 Soil salinity prediction by the MLR models 
3.4.1 Prediction models for ECe 


The prediction models for ECe and the fitting parameters used in the MLR models and their 
respective R? are presented in Table 2. All MLR models had high R? values ranging from 0.78 to 
0.91 for the surface soil layer and from 0.77 to 0.95 for other soil layers (P<0.05). For the 
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ECa-gm31 data, the lowest R? of 0.77 was calculated for the 1.2-3.0 m soil layer, and the highest 
R? was observed for the 0.0-0.6 m surface soil layer. 


Table 2 Multiple linear regression coefficients of predicted electrical conductivity of saturated soil-paste extract 
(ECe) at different soil depths 


p 


EM device Soil depth (m) a b c R N 
0.0-0.6 0.051 0.005 —6.37 0.91 
0.6-1.2 —0.003 0.017 17.92 0.78 

EM31 6 
1.2-3.0 —0.050 0.100 35.80 0.77 
3.0-6.0 0.048 0.003 19.63 0.95 
0.0-0.6 0.006 0.012 3.20 0.86 

EM38 23 
0.6-1.2 0.003 0.038 6.21 0.85 


Note: The results were obtained considering that ECa-gmy and ECa-emu were calibrated for the average ECe measured at 23 sampling sites 
with EM38 (0.0-1.2 m soil layer) and 6 sampling points with EM31 (0.0-6.0 m soil layer). N means the total number of sampling sites. 


3.4.2 Prediction models for ECe 

Prediction of ECe using the MLR models showed an increase in soil salinity with depth (Table 3). 
From spring to autumn 2015, the highest mean ECe was found in the soil layer that was 1.2-3.0 m 
deep with a maximum value of about 59.90 dS/m in spring; the lowest salinity of 10.02 dS/m was 
observed in the surface soil layer at 0.0-0.6 m in autumn 2015. During autumn 2015, the 
maximum of ECe was 52.71 dS/m in the soil layer with a depth of 1.2—3.0 m (mean value 42.82 
dS/m), while in autumn 2016, the highest value of 55.39 dS/m was estimated for the soil layer 
with a depth of 1.2-3.0 m. The highest mean ECe (46.44 dS/m) was predicted for the deeper soil 
layer at 3.0—6.0 m. 


Table 3 Descriptive statistics of predicted ECe and the estimation of ECe using Ordinary Kriging method at 
different soil depths in spring 2015, autumn 2015, and autumn 2016 


Soil Predicted ECe Kriging estimation of ECe 
Period depth (m) Mean SD Max. Min. Median kurtosis Skew Mean SD 
(dS/m) (dS/m) (dS/m) (dS/m) (dS/m) (dS/m) (dS/m) 
0.0-0.6 13.47 1:52 18.02 10.50 13.27 —0.38 0.31 14.02 1.17 
0.6-1.2 23.31 2.57 31.87 18.40 23.06 0.10 0.48 24.98 2.08 
Spring 2015 
1.2-3.0 43.26 3.49 59.90 32.36 41.58 2.52 1.14 43.02 1.77 
3.0-6.0 42.32 1.70 46.69 38.63 42.87 —0.25 0.31 42.30 0.83 
0.0-0.6 12.74 1.27 18.29 10.02 12.53 0.10 —0.29 12.80 1.01 
0.6-1.2 19.72 1:23 23.39 16.89 20.02 —0.09 —0.09 22.55 2.27 
Autumn 2015 
1.2-3.0 42.82 2.54 52.71 36.13 41.45 0.75 0.62 42.52 1.43 
3.0-6.0 42.21 1.68 46.78 38.61 43.03 —0.24 0.38 42.10 0.88 
0.0-0.6 13.74 1.32 17.79 11.25 13.74 0.47 0.95 12.60 0.56 
0.6-1.2 23.11 2.60 32.39 19.38 22.87 1.48 1.16 21.50 1.62 


Autumn 2016 
1.2-3.0 44.36 3.02 55.39 37.96 44.23 0.93 0.56 44.74 1.15 


3.0-6.0 46.44 1.12 49.31 43.03 46.59 0.18 —0.00 46.46 0.50 


3.5 Spatial and temporal variability in soil salinity 

3.5.1 Variogram analysis 

All the variograms show a positive nugget of ECe ranging from 0.70 to 14.72 dS/m (Table 4). The 
ratio of variance and structural variance related to the sill, MAE, and RMSE are presented in 
Table 4. These results indicated a 39.30% spatial variation in EC in the surface soil layer and a 
60.00% variation in the 6.0-m-deep soil layer in spring 2015. However, this variability doubled in 
autumn 2016 in the surface soil layer (0.0—0.6 m). 
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Table 4 Variogram characteristics and cross-validation statistics for ECe at different soil depths in spring 2015, 
autumn 2015, and autumn 2016 


Variogram of ECe Cross validation of ECe 
Period Soil depth (m) 
Nugget (dS/m) Sill (dS/m) Range (m) Nugget/Sill MAE (dS/m) RMSE (dS/m) 

0.0-0.6 1.12 2.85 27.59 0.39 —0.006 1.22 

0.6-1.2 3.23 10.67 48.00 0.30 —0.006 2.33 
Spring 2015 

1.2-3.0 14.72 23.55 26.27 0.63 0.010 4.29 

3.0-6.0 2.46 4.09 12.14 0.60 0.006 1:57 

0.0-0.6 1.14 3.00 33.76 0.38 —0.005 1.39 

0.6-1.2 5.90 14.20 37.95 0.42 —0.002 2.78 
Autumn 2015 

1.2-3.0 7.02 10.48 45.25 0.67 0.002 2.69 

3.0-6.0 2.06 3.94 16.00 0.52 —0.007 1:57 

0.0-0.6 1.60 2.34 48.00 0.68 0.000 1.39 

0.6-1.2 5.70 12.53 96.00 0.45 —0.007 3.04 
Autumn 2016 

1.2-3.0 5.43 7.75 64.00 0.68 —0.009 2.49 

3.0-6.0 0.70 1.22 15.67 0.57 0.001 0.98 


Note: MAE, mean average error; RMSE, root mean squared error. 


The exponential semi-variogram was found to be the best method to estimate ECe because it 
can predict soil ECe values for non-sampled sites. Figure 4 shows the exponential 
semi-variograms of the measured data and the fitted MLR models. The results demonstrated that 
the efficiency ranged from 12.10 to 48.00 m in spring 2015 and reached 96.00 m in autumn 2016, 
which indicated that the smallest distance from EC sampling site was 12.1 and 48.0 m in spring 
2015 and ranged between 15.8 and 96.0 m in autumn 2016. 

3.5.2 Prediction accuracy of the Ordinary Kriging method 

Performance indicators for ECe mapping for the different soil layers using the Ordinary Kriging 
method are given in Table 4. The RMSE ranged from 0.98 to 4.29 dS/m. The lowest RMSE was 
noted for the last sampling in autumn 2016 in the deep soil layer (3.0-6.0 m), while the highest 
RMSE value was observed for the 1.2-3.0 m soil layer sampled in spring 2015. A high RMSE 
(>2.00 dS/m) was observed in both layers between 0.6 and 3.0 m in all three sampling periods. 
However, the MAE was extremely low and did not exceed 0.010 dS/m (Table 4). The 
cross-validation statistics of inverted soil ECe in the soil layers of 0.0-0.6 and 0.0-1.2 m are given 
in Table 5. The RMSE ranged from 1.73 to 3.67 dS/m. The lowest RMSE was noted for the 0.0-1.2 
m soil layer, while the highest RMSE value was observed for the 0.0-0.6 m surface soil layer. 


3.6 Spatial and seasonal maps of soil salinity 


Maps of predicted ECe are shown in Figure 5. The ECe estimates by the Ordinary Kriging 
method ranged between 8.00 and 16.00 dS/m in the 0.0-0.6 m deep surface soil layer and ranged 
between 16.00 and 32.00 dS/m in the soil layer between 0.6 and 1.2 m. The estimates reached a 
maximum of 52.00 dS/m in the deep soil layer (3.0-6.0 m). The inverted ECe maps are shown in 
Figure 6, where the inverted ECe varied between 8.00 and 24.00 dS/m in the 0.0-0.6 m surface 
soil layer and between 16.00 and 22.0 dS/m in the 0.0-1.2 m soil layer. An equivalent difference 
was observed among the three sampling periods. The ECe values for the surface soil layer 
decreased from spring 2015 to autumn 2016 but increased for the deep soil layer (Figs. 5 and 6). 

To quantify the temporal variability in soil salinity, CV values for each ECa-gm measured over 
time between the first and second sampling periods from both conductivity meters are shown in 
Figure 7. These maps show that soil salinity varied from one soil depth to another, depending on 
the conductivity meter used. The ECa-emssų data had the highest CV. The variation in ECa-em3gu 
at a depth of 0.6 m and 37.00% of zoning was greater than that with 25.00% of zoning, especially 
at the level of the northwest of the study area. This variation decreased to 5.00% at a depth of 1.2 
m (ECa-emssv; Fig. 7). 
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Fig. 4 ECe variograms predicted by the multiple linear regression (MLR) models in spring 2015 (a), autumn 


2015 (b), and autumn 2016 (c) 


Table 5  Cross-validation statistics of ECe inverted at different soil depths in spring 2015, autumn 2015, and 


autumn 2016 


Cross-validation statistic 


Period Soil depth (m) 
MAE (dS/m) RMSE (dS/m) 

0.0-0.6 2.10 2.72 
Spring 2015 

0.0-1.2 1.37 1.73 

0.0-0.6 1.81 2.34 
Autumn 2015 

0.0-1.2 1.66 2.28 

0.0-0.6 2.66 3.67 
Autumn 2015 

0.0-1.2 1.65 2.06 


The variation in ECa-em3ın (4.0 m) remained moderately stable between 10.00%—15.00% for 
45.00% of zoning for the northeast of the study area; a similar observation was obtained for 
ECa-emaiv (6.0 m) for only 38.00% of zoning. The skewness for all measurements was between 
0.52 and 1.11, and presented a symmetrical distribution. However, the skewness for ECa-Em3sv 
was higher than 1.00, so the distribution is asymmetric (Table 6). 
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Fig. 5 Spatial and temporal distributions of the one dimensional inversion in ECe for the soil layers of 0.0-0.6 
m (a, b, c) and 0.0-1.2 m (d, e, f) in spring 2015, autumn 2015, and autumn 2016 


4 Discussion 


4.1 Soil salinity monitoring over space and time 


Temporal variations of soil salinity were mainly found in the southeast and northeast of the study 
area. High variability was particularly evident in the surface soil layer, especially during the 
second sampling period conducted in autumn 2015. This is attributed to the heavy rainfalls prior 
to sampling. The temporal dynamics could have been determined by drip irrigation practices or 
rainfall fluctuations during the winter and summer seasons. Heavy rainfall events occurred on 16 
October 2016 (>50.0 mm) and 11-12 July 2015 (about 100.0 mm) that may have caused 
substantial salt leaching. 

There was a distinctive difference in the predicted ECe values between the northwest axis in 
the center and the eastern and western boundaries. The maps of estimated ECe showed a great 
difference in ECe of the surface and deep soil layers, where ECe values increased with depth in 
all sampling periods, as mentioned earlier (Zhang et al., 2021). The ECe decreased in the surface 
soil layer between the first and second sampling periods and increased during the last sampling 
period (autumn 2016). This decrease is likely due to the clayey texture of soils and the input of 
water from rainfall and irrigation, which helped to maintain a high GWC in the surface soil layer. 
In autumn 2016, capillary action toward the surface and the accumulation of salts at the surface 
horizons led to surface salinization, as reported previously (Hachicha et al., 2000; Doolittle and 
Brevik, 2014). 


4.2 Model performance and soil salinity mapping 


The methods developed by Lesch et al. (1998) to predict ECe from EM38 measurements using 
MLR models were applied to a single field (Lesch et al., 1998; Amezketa, 2007). Spatial analysis 
is generally complicated due to the following reasons: (1) the inability to acquire deep-core 
samples from dry, hard, and compact clay soils; (2) the extremely high spatial and temporal 
variability in the soil GWC that influences ECa and the soil salinity variability in space and time; 
and (3) the seasonal variability in the salinity of surface soil layer (Safar, 1983). Our salinity 
mapping approach was based on coupling the standardized laboratory method of ECe 
determination with ECa measurement using geophysical methods (EM31 and EM38) (Zhao et al., 
2020). 


202203.00037v1 


chinaXiv 


ChinaXivA ERAT! 


JOURNAL OF ARID LAND 


Spring 2015 Autumn 2015 Autumn 2016 


(al) 0.0-0.6 m (a2) 0.0-0.6 m (a3) 0.0-0.6 m 


Legend 
ECe (dS/m) 
E <8 
E 8-12 
E 12-14 
EE 14-16 
E 16-18 
E 18-20 
= 20-22 
E 22-24 
m 24-28 
= 28-32 
32-36 
E 3640 
E 40-42 
E 42-44 
Mi 44-46 
Mm 46-48 
Wl 48-52 
— Isoline 
(cl) 1.2-3.0m (c2) 1.2-3.0 m (c3) 1.2-3.0 m % Measurement 
point 
0 30m 
bl 


(b1) 0.6-1.2 m (b2) 0.6-1.2 m (b3) 0.6-1.2 m 


(d1) 3.0-6.0 m (d2) 3.0-6.0 m (d3) 3.0-6.0 m 


Fig. 6 Spatial and temporal distributions of ECe obtained from measurements of ECa-gm3s in the soil layers of 
0.0-0.6 m (al—a3) and 0.6-1.2 m (b1—b3) and ECa-rmsi in the soil layers of 1.2-3.0 m (cl—c3) and 3.0-6.0 m 
(d1—d3) m in spring 2015, autumn 2015, and autumn 2016 


The steps of the methodology we adopted are as follows: first, the analysis with the EM38 and 
EM31 measurements could be extended over a large area and to deeper soil layers of. Second, the 
predicted salinity at different soil depths was used to analyze the vertical salinity profile (Michot 
et al., 2013). The MLR models were applied to estimate the variability in EC and describe the 
spatial variability in soil salinity with reasonable accuracy (1.20 dS/m<RMSE<4.30 dS/m). 
Findings showed that over a period of three years, the spatial pattern of soil salinity had not 
changed much, especially in the deep soil layer (3.0—6.0 m). 
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Fig. 7 Coefficient of variation (CV) of ECa-emn (a and b) and ECa-emv (c and d) measurements using EM31 
and EM38 devices 


Table 6 Descriptive statistics of the CV of ECa in spring 2015, autumn 2015, and autumn 2016 


CV of ECa 
Measurement profile Device 
Mean (%) Median (%) SD (%) Kurtosis Skewness 

Vertical 10.20 9.07 6.20 1.18 1.11 
EM38 

Horizontal 21.60 21.21 8.81 —0.84 0.12 

Vertical 11.20 11.11 3.93 —0.36 0.14 
EM31 

Horizontal 12.80 12.37 5.62 —0.03 0.52 


The need to adapt the regression models to local soil context was indicated by the lower 
precision of the general linear regression model (which integrated the ECa-emu and ECa-emv 
measurements) as compared to the MLR models, despite the limited number of sites used in the 
latter. As described by Triantafilis et al. (2001b), EM38 and EM31 measurements can be used for 
a more accurate prediction of ECe, as they allow more samples to be measured cost-effectively 
and without disturbing the soil structure. Cross-validation can provide evidence of the efficiency 
of these conductivity meters to predict ECe at the plot scale as per the findings of Corwin and 
Lesch (2005), Omonode and Vyn (2006), and Huang et al. (2016). 


5 Conclusions 


This work aimed to link two conductivity meters for the measurements of EC to assess the spatial 
and temporal variations in soil salinity in an irrigated field. Variation in soil salinity of alluvial 
clay soil was determined in terms of ECe measurements performed using two conductivity 
meters, i.e., EM38 and EM31. An MLR model built using the measured soil ECe and ECa-gm3 
was found to be a good predictor of ECe for the surface soil layer and an MLR model built using 
ECa-gm31 was effective to predict ECe for the deep soil layer (about 6.0 m). 
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The ECe maps obtained using the Ordinary Kriging method showed that spatial changes in 
salinity remain relatively constant; ECe dropped in the surface soil layer from spring 2015 to 
autumn 2016 but increased in the deep soil layer. All maps showed a low CV of ECe (<10.00%) 
in the center of the study area. The variation in salinity at the soil surface (ECa-zm3gn at the depth 
of 0.6 m) was greater than that at the depth of 1.2 m (ECa-gm3sv at the depth of 1.2 m). At the soil 
surface, the variation in salinity was greater than 25.00%, while at a depth of 1.2 m, it was about 
5.00%. The variation in ECa-gm3in (4.0 m) remained moderately stable between 10.00% and 
15.00% for 45.00% of zoning in the northeast of the study area; similar variation was recorded in 
ECa-emaiv (6.0 m) for only 38.00% of the zoning area. These methods helped to track salinity 
trends and their significance in the field; they can prove useful for the identification of spatial and 
temporal features, leading to more effective management of agricultural soil, regardless of tillage 
or irrigation conditions. From our standpoint, the application of EM38 and EM31 was effective in 
Central Tunisia for large-scale monitoring of salinity in the soils. 
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